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Abstract 



We present a new parton level Monte Carlo program for the calcula- 
tion of jet cross sections in Deep Inelastic Scattering based on Born 
and next-to-leading order matrix elements. Using a class of invari- 
ant jet definition schemes, the program allows for the calculation of 
differential distributions of jet cross sections in the basic kinematical 
variables (like s, x, y, W 2 , Q 2 . . .) as well as for total jet cross sections. 
Various kinematical cuts can be chosen from an input file. 



1 Introduction 



The start-up of the HERA electron-positron collider marked the beginning of a 
new aera of experiments exploring Deep Inelastic Scattering (DIS) of electrons and 
protons. One of the topics to be studied at HERA will be the Deep Inelastic (= 
high Q 2 ) production of multi jet events, where good event statistics are expected 
allowing for precision tests of QCD. 

In this paper, we present a Monte Carlo program DISJET that allows for 
the calculation of next-to-leading order (NLO) (1 + 1),(2 + 1) and leading order 
(LO) (3 + 1) jet cross sections^] in DIS. Using a class of invariant jet definition 
scheme [0 the program allows for the calculation of differential distributions of 
jet cross sections in the basic kinematical variables as well as for the calculation 
of total jet cross sections. 

The source code is written in Fortran and running ability is tested on a 
DEC-Station and an ALPHA-Station running the Operating Systems VMS and 
on a Sun-Sparc station running under UNIX. The Monte Carlo routines are using 
the VEGAS-package for numerical integration of functions depending on up to 
ten variables. Parton distributions are incorporated from the PDFLIB-package 
H. Each (n + 1) jet cross section in DISJET is calculated in fixed order QCD 
perturbation theory. Therefore a comparison of jet-measurements with our pre- 
dicted QCD-results provides a direct tool of determining a s or Ams. 

Changes to DISJET 1.0 

An error in the finite part of the 0(a 2 ) virtual corrections to the quark initiated 
subprocess is corrected. An error in the labeling of the 0(a 2 ) (3+l)-jet subpro- 
cesses [iproc = 1001,1002,1003] is corrected. An additional cut in the variable 
z = — , where p\ is the four momentum of one of the final partons in (2+1) jet 
production is implemented. 

2 Matrix Elements 

In DISJET complete fixed order QCD matrix elements for the process 

e - (Z) + proton(P) — > proton remnant(p r ) + parton l(pi) . . . + parton n(p n ) (1) 

(n = 1, 2, 3) are used. Reaction (H) proceeds via the exchange of an intermediate 
vector boson V = 7*, Z, W. At present only the exchange by a virtual photon 
is incorporated in DISJET. We denote the 7*-momentum by q, the absolute 
square by Q 2 , the center of mass energy by s, the square of the final hadronic 

1 " + 1" denotes the remnant jet 
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mass by W 2 and introduce the scaling variables x and y: 



q = I -I' 

Q 2 = -q 2 = xys > 

s = {P + lf 

W 2 = P) = (P + qf (2) 

o 2 

x = — (0 < x < 1) 
2Pq K ~ ' 

Pq 

y = -f x (0<y<l) 

At fixed s, only two variables in (^) are independent, since e. g. 

X W 2 = (1 - x)Q 2 , Q 2 = xys. 

Averaging over the azimuthal angle between the parton plane and the lepton 
plane (I, I') (in the (7*-initial parton)-cms) the jet cross sections for the exchanged 
virtual photon factorize the following y dependence [0: 

der[n-jet] = {(1 + (1 - y) 2 ) da u+L [n - jet] - y 2 da L [n - jet]} (3) 

In eq. (|3|) <Ju+L[ n — jet] and a^n — jet] denotes the cross section contributions 
from an unpolarized and longitudinal polarized gauge boson. These helicity cross 
sections are linearly related to polarization density matrix elements of the virtual 
boson. One has: 

a u = + h — ) cTL = h 00 a u+L = au + a L (4) 

where h mm i = e*^(m)H fJ,u e u (m / ), (m,m f = +,0, — ) and e M (±)(e M (0)) are the trans- 
verse (longitudinal) polarization vectors of the 7* in the (7*— initial parton)-cms. 
The helicity cross sections can technically be obtained by the following projections 
on the (partonic) hadrontensor H^, which is calculated in fixed order perturba- 
tion theory, 

/ 1 3x \ 
ov + l = (-3^ + — PmM # M >-jet] (5) 

ol = -^v^VvE^\n- jet] (6) 
pq 

where p denotes the momentum of the initial state parton and x p = Q 2 /(2pq). 
The available jet multiplicities in DISJET are listed below: 

1) (1+1) jet 0(a° 8 ) and 0(a]) 
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2) (2+1) jet 0(al) 

3) (2+1) jet 0{a 2 ) 

4) (3+1) jet 0(a 2 ) 

The physics of 1) and 2) is extensively discussed in Q. For numerical results 
see also [|| Matrix elements of the complete contributions to (2+1) jet 0(a 2 ) 
corrections are first discussed in ||. In DISJET, a modified version of these 
matrix elements is implemented including the full NLO scale dependence. NLO 
numerical results for (2+1) jet cross sections and rates (based on DISJET) in- 
cluding the contributions from all helicity cross sections are presented in |J . The 
(2+1) jet NLO contributions originating from the projection with the metrical 
tensor — (see eq. (|5])) on the hadron tensor and a discussion of the results are 
given in @. Results for the 0(a 2 ) (3+1) jet cross sections are discussed in detail 
in §. 

The general structure of the hadronic jet cross sections within the framework 
of perturbative QCD is given by 

da had [(n+l)-jet} = J dq f a (rj, M 2 ) da a (p = V P, a s (M 2 )) (7) 

where one sums over a = q,q,g. f a {r],M 2 ) is the probability density to find a 
parton a with fraction r\ in the proton if the proton is probed at a scale M 2 . a a 
denotes the partonic cross section from which collinear initial state singularities 
have been factorized out (in the NLO calculation) at a scale M 2 . Note that a a 
depends in NLO also explicitly on the renormalization scale M r and factorization 
scale Mf. DISJET allows the user to choose the mass scales M r and Mf in eq. 
(0) according to the following input variables a x , b x and c x \ 

M 2 = a x Q 2 + b x W 2 + c xP 2 T (8) 

where x = f,r (see also next section). In eq. (|j) pr denotes the transverse 
momentum of the two partonic jets in (2 + 1) jet production. The transverse 
momentum is defined with respect to the 7* direction. c x 7^ can only be chosen 
in the case of (2+1) jet production, otherwise c x is set to zero. Furthermore, to 
avoid to small scales for perturbation theory, the renormalization and factoriza- 
tion scales are cut at 2 GeV 2 , i.e. 

M 2 = max{2 GeV 2 , M 2 } (9) 

In order to calculate the (n + 1) jet cross sections we have to define what we 
call (n + 1) jets by introducing a resolution criterion. As has been elaborated 
in detail in |], || energy angle cuts are not suitable for an asymmetric machine 
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with its strong boosts from the hadronic cms to the laboratory frame. As a jet 
resolution criterion we use the invariant mass cut criterion introduced in Jl], 
such that 

Sij = (Pi + Pj) 2 > M 2 = max {y cut M 2 , M 2 } (i,j = l,...,n,r; i # j) 

(10) 

where y cut is the resolution parameter which should be chosen in [0.01 . . . 0.08] 
(see fig. 5 in M c is a typical mass scale of the process which defines the jet 
definition scheme. This mass scale can be chosen by the user according to the 
values of a c , b c and a constant C [GeV 2 ]: (see also section 3) 

M 2 C = a c Q 2 + b c W 2 + C (11) 

The extreme choices a c = 0(b c = 0) and C = correspond to the "W — 
scheme" ("Q — scheme") as introduced in [p]]. M in eq. ( |10|) is an additional fixed 
mass cut which we have introduced in order to clearly separate the perturbative 
and non-perturbative regime in the case, where M 2 is small (i.e. Q 2 and/or W 2 
and C small). A reasonable choice for Mq is Mq = 2 GeV 0. 



3 Input parameters 

The input parameters for DISJET are written to a file referred to as unit 4 by 
the local Fortran compiler, i.e. for004.dat on VMS and fort. 4 on UNIX. 

Parameters to choose are: (an example for the input file is given in the ap- 
pendix) 

• srs: real, f6.2, [GeV] 

center of mass energy yj \P + I) 2 

• ias: integer, il 

choose one-loop ore two-loop formula for a-j^g- The value of a s is matched 
at the thresholds q = m q and the number of flavours rif in a s is fixed by 
the number of flavours for which the masses are less than M r . Furthermore 
the numbers of quark flavours that can be pair-produced are set equal to 
n/ chosen in a s . 

1: 1-loop- formula 
2: 2-loop- formula 

• ilambda: {0, 1} integer, il 

0: in a s is chosen consistent to the value in the chosen parton dis- 

tribution functions (see below). 
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1: Ajlg in a s can be chosen by hand (—dlam), see next item. 

- dlam: real, f5.4, [GeV] 

only relevant in the case ilambda = 1, otherwise this number is changed 
in the main program. 

If ilambda = 1 , Aj^=dlam, real, f5.4., [GeV]. 

• iaem: {0, 1} integer, il 

fixed or Q 2 dependent a e i ectromagnetic . 

0: a = 0.00729735 fixed 
1: Q 2 dependent a 

• icross: {0, 1} integer, il 

0: calculate the jet cross section as defined in (|3|) in [pb] 

1: only effective for iproc < 1000 (see below) 

calculates the helicity cross sections <Ju+l jet] or &L[{n+l) — jet] 
for fixed x and Q 2 (n = 1,2) 

- x: [0.0001,1]: real, f5.4 

If icross=0, x is changed by the integration routine. 

- Q"2: [AGeV 2 ,s]: real, fl0.2 
same as x 

• ipola: {10, 11, 12}: integer, i2 

choose helicity cross sections in the case of (n+1) jet production (n = 1,2). 
Depending on icross, the full cross sections (including all coefficients) or 
the helicity cross sections cr u+LjL (for fixed x, Q 2 ) are calculated. 

10: take sum of all helicity cross sections, i.e. calculates da[n — jet] as 
given in (|3[) (only possible for icross = 0). 

11: calculates only the term ~ dau + i,[n — jet] in ([$]) 

12: calculates only the term ~ da^n — jet] in (||) 

• iproc: {10, 100, 101, 102, 190, 191, 192, 1001, 1002, 1003}: integer, i4 
characterization of the process 

10: (1+1) jet = a tot in O(a s ) 

100: (2+1) jet: complete, i.e. sum of quark and gluon initiated processes 

101: (2+1) jet: quark initiated 

102: (2+1) jet: gluon initiated 
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190: 0(a^)-contribution to o tot in the DIS-scheme: complete 

191: 0(a*)-contribution to o t ot in the DIS-scheme: quark initiated 

192: 0(a*)-contribution to a to t in the DIS-scheme: gluon initiated 

1001: (3+1) jet: q -> qGG (class D (§) 

1002: (3+1) jet: G -> Gqq (class C [|) 

1003: (3+1) jet q -»■ ggY (class E J§) 

» istruc: {0, . . . , 7}: integer, il 
only relevant if iproc= 100, 101, 102 

0: complete NLO corrections, (note: ias should be set equal to 2) 

1: Leading order (0(a s )) 

2: NLO virtual corrections (finite part) 

3: NLO final state real corrections 

4-7: NLO initial state real corrections: 

4: delta-part 5: hard-part 6: plus-part 7: substracted-part 
(note: ias should be set = 2 for istruc=2-7) 

» y_c: [0.01,0.1]: real, f5.4 
jet resolution parameter y cut as defined in eq. (]T0|) 

» M.(T2r. real, f3.1, [GeV 2 ] 
additional fixed jet cut mass as defined in eq. ( [10]) 

» M_f 2 = cif Q" 2 + bf W" 2 + Cf pf,: the factorization mass in eq. (0) can 
be chosen by the three values a/,&/ and c/ in a format f4.2 (c/ ^ only 
effective for iproc = 100, 101, 102) 

» M_r" 2 = a r Q" 2 + b r W" 2 + c r p\.\ the renormalization mass in eq. (0) 
can be chosen by the three values a r , b r and c r in a format f4.2 (c r 7^ only 
effective for iproc = 100, 101, 102) 

» M_c 2 = a c Q" 2 + b c W 2 + C [GeV] 2 : the jet cut mass as defined in eq 
(|10"D can be chosen by the three values a c , b c and C in a format f4.2 



2 rp mi -. real, f6.2 



in in 



Pt 

minimal required transverse momentum for the jets (transverse with respect 
to the 7* direction). This entry is effective only for iproc = 100, 101, 102. 
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[0.,1.]: real, f8.2 

Cut on the invariant z = —, where pi is the four momentum of one of 

pq ' rL 

the final partons in (2+1) jet production. This entry is effective only for 
iproc = 100, 101, 102. 

• Nptype, Ngroup, Nset integer, i3,i3,i3 

choose parton distribution functions according to 0. 

• iacc: {1, 2, 3} integer, il 

some predefined settings for the accuracy the VEGAS-integration: 
1: acc = 0.05, ncall = 20000 
2: acc = 0.005, ncall = 200000 
3: acc = 0.0005, ncall = 1000000 

• x m i n ,x max : [0.0001,1] real, f8.6, Integration limits 

• VminiVmax- [0,1] real, f8.6, Integration limits 

• QmimQmax- [ 4 > s ] real > f8 - 2 > [GeV 2 ] , Integration limits 

• ^min^max- [5, s] real, f8.2, [GeV 2 ], Integration limits 

• The next integer (format i2) determines the number of single differential 
cross sections to be calculated. 

• the following input values are data-cards for the calculation of single dif- 
ferential cross sections (see the last column). The user may change the 
entries in the first two columns. The first (second) column defines the lower 
(upper) kinematical limit of the variable specified in the last column. The 
user may also change the single differential cross sections to be calculated 
by suitable changings in the subroutine dsigma. 



4 Output 

The output of DISJET is written to the standard output as well as to a file 
unit 98. The output file contains all input values specified in the unit 4 input file 
and all results calculated by VEGAS, i.e. the integrated cross section and the 
differential cross sections. 

The differential cross sections will be written to file unit 98 in the following 
format. 

• the (up to) 10 differential cross sections are written to the file one after 
another 
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• every package of data belonging to one differential cross section is preceded 
by the string given in the data cards of the input file 

• each differential cross section is presented line by line for each requested bin 
in the form: 

mean value of the variable in the bin, differential cross section in this bin, 
number of points in this bin. 

5 Summary and conclusion 

We present the first version of the Monte Carlo program DIS JET. The program 
allows to calculate NLO (1+1) and (2+1) jet cross sections and LO (3+1) jet cross 
sections in DIS. The contributions from all helicity cross sections are included. 
Note that DIS JET is not an event generator like LEPTO or HERMES [§. 
Furthermore, the freedom to choose a jet definition scheme is restricted to the 
schemes discussed in the text and does not include the jet definition proposed in 

!- 

Some improvements and extensions are planed for further versions of the pro- 
gram. We would appreciate comments and suggestions on the program to be send 
to mirkes@phenom.physics.wisc.edu. 

Acknowledgements: 

It is a pleasure to thank I.H. Park and J. Rathsman for their efforts to test the 
MC program DISJET and helpful suggestions concerning some extensions of 
the program. This work is supported in part by the U.S. Department of En- 
ergy under contract No. DE-AC02-76ER00881, and in part by the University 
of Wisconsin Research Committee with funds granted by the Wisconsin Alumni 
Research Foundation. 

References 

[1] J.G. Kdrner, E. Mirkes, G. Schuler: Int. J. Mod. Phys. A4 1989 1781. 
[2] G. P. Lepage, J. Comp. Phys. 27 (1978) 192. 
[3] H. Plothow-Besch, CERN-PPE/ 92-123 

[4] T. Brodkorb, J.G. K6rner, E. Mirkes, G. Schuler: Z. Phys. C 44 (1989) 
415-432. 

[5] T. Brodkorb and E. Mirkes, Madison preprint, MAD/PH/820. 
[6] T. Brodkorb, J.G. Kdrner: Z. Phys. C 54 (1992) 519. 



9 



[7] D. Graudenz Phys. Lett. B 256 (1992) 518 
LBL-34147 preprint, July 1993; 

[8] Proceedings of the HERA Workshop 1991, Vol. 1-3, DESY (eds. W. 
Buchmuller and G. Ingelman.) 

[9] S. Catani, Yu.L. Dokshitzer and B. Webber, Phys. Lett. B 285 (1992) 291 
S. Catani, Yu.L. Dokshitzer and B. Webber, Phys. Lett. B 322 (1994) 263 



10 



Appendix 



Example for an input-file 
for004.dat under VMS 
fort. 4 under UNIX 



SRS= 295. 
IAS= 1 

ILAMBDA= 0, DLAM= .1970 
IAEM= 1 

ICR0SS= 0, X= .010, Q~2= 10000.00 
IP0LA= 10 
IPR0C= 100 
ISTRUC= 1 

Y_C= .0200, M_0~2= 4.0 GEV~2 

M_F~2= 1.00 Q~2 + 0.00 W~2 + 0.00 PT"2 
M_R~2= 1.00 Q~2 + 0.00 W~2 + 0.00 PT~2 
M_C~2= 0.00 Q~2 + 1.00 W~2 +0.00 
PT2MIN= . 

NPTYPE = 1, NGROUP = 3, NSET =31 
IACC= 1 

XMIN= 0.001000 XMAX= 0.100000 
YMIN= 0.040000 YMAX= 0.950000 
ZP.MIN = 0.00 ZP.MAX = 1.00 
Q~2_MIN= 10.00 Q~2_MAX= 87025.00 
W~2_MIN= 600.00 W~2_MAX= 87025.00 
10 

-3., -0.1, 32,0,0, 'LG(X) ' 

3.D-4, l.D-1, 32,0,0, 'X' 

2., 10. DO, 32,0,0, J SQRT(Q"2) ' 

10., 1.D5, 32,0,0, 'Q~2' 

3., 200. DO, 32,0,0, J SQRT(W"2) ' 

10., 1.D5, 32,0,0, 'W~2' 

-2., -0.1, 32,0,0, J LG(Y) ' 

0.05, 0.95, 32,0,0, 'Y' 

0.67, 0.99, 32,0,0, 'THRUST' 

0. , 50., 32,0,0, 'PT' 
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